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ABSTRACT 


Thermal  instabilities  in  a  viscoelastic  rod  under  cyclic  loading 
are  discussed  by  determining  the  stresses  and  temperature  in  a  visco¬ 
elastic  rod  insulated  on  its  lateral  surface  and  driven  by  a  sinusoidal 
stress  at  one  end.  Temperature  dependence  of  the  complex  Young's 
modulus  of  the  rod  and  the  effect  of  thermomechanical  coupling  are  in¬ 
cluded  in  the  analysis.  A  method  of  finite  differences  is  used  to 
directly  determine  the  steady  state  stresses  and  temperature  without 
obtaining  the  complete  time  history  of  the  process.  The  iterative 
algorithm  used  is  very  useful  and  converges  rapidly  for  a  wide  range 
of  driving  stress  amplitudes  and  frequencies.  It  is  found  that  rapid 
rise  of  temperature  to  dangerous  levels  occurs  for  relatively  low 
values  of  driving  stress  amplitudes,  especially  if  the  driving  fre¬ 
quency  is  close  to  one  of  the  critical  frequencies  of  the  rod.  Dras¬ 
tic  softening  of  the  rod  leads  to  large  strains.  Thus,  failure  of  the 
rod  could  occur  at  low  values  of  the  driving  stress. 
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INTRODUCTION 

Thermal  effects  In  viscoelastic  materials  under  cyclic  loading  are 
of  great  practical  importance.  These  materials  are  dissipative  in  na¬ 
ture.  Continuous  dissipation  of  mechanical  energy  into  heat  coupled 
with  the  fact  that  they  are  poor  conductors  of  heat  can  lead  to  very 
high  temperatures  inside  the  materials  and  consequent  failure.  In  ad¬ 
dition,  the  mechanical  properties  of  the  materials  are  strongly  temper¬ 
ature  dependent,  the  complex  elastic  modulii,  in  general,  being  in¬ 
versely  proportional  to  some  power  of  the  temperature.  This  rapid 
softening  of  the  material  often  leads  to  large  strains  for  relatively 
small  values  of  stress. 

Several  authors  have  studied  thermal  effects  in  viscoelastic  media. 

"fc 

Tormey  and  Britton  [1]  carried  out  vibration  tests  of  solid  propellant 
rocket  motors.  Heating  due  to  vibration  near  a  resonant  frequency  for 
several  hours  caused  the  material  to  soften  to  such  an  extent  that  some 
of  it  flowed  out  of  the  motor.  Hunter  [3]  derived  a  set  of  thermomechan- 
ically  coupled  equations  for  the  propagation  of  stress,  strain  and  tem¬ 
perature  fields  in  viscoelastic  solids.  He  assumed  thermorheologically 
simple  behavior  and  used  a  double  time  integral  expression  for  the  dis¬ 
sipation  function.  Petrof  and  Gratch  [3]  considered  longitudinal  wave 
propagation  in  a  finite  rod.  They  assumed  that  the  material  behaves 
like  ideal  rubber  so  that  deformation  occurs  without  any  change  in  po¬ 
tential  energy.  They  used  an  integral  form  of  the  constitutive  equa¬ 
tions.  The  cransient  response  of  the  same  problem  was  studied  by 

•jf 

Numbers  in  square  brackets  designate  references  at  the  end  of  this 
report. 
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Wolosewick  and  Gratch  [4].  They  found  that  due  to  the  high  damping  in 
viscoelastic  materials,  the  mechanical  transients  decay  very  quickly, 
and  a  steady  state  of  oscillation  with  slowly  increasing  temperature  can 
be  reached  within  one  or  two  cycles. 

In  a  series  of  papers  [5,6,7],  Schapery  derived  field  equations  for 
viscoelastic  media  with  thennomechanical  coupling  and  solved  some  prob¬ 
lems.  In  [6]  he  studied  the  problem  of  steady  periodic  shear  oscilla¬ 
tions  using  the  complex  modulus  form  of  the  constitutive  equations.  He, 
however,  neglected  inertia  and  his  application  was  limited  to  lower  fre¬ 
quencies  than  considered  here.  In  [7]  he  developed  and  used  variational 
principles  to  solve  problems  involving  bodies  that  were  either  massless 
or  with  concentrated  mass,  and  in  his  last  example  of  a  ‘solid  cylinder 
with  distributed  mass1  he  only  gives  a  first  approximation  to  the  solu¬ 
tion.  This  paper  1.7]  has  been  discussed  in  more  detail  in  [25]. 

Huang  and  Lee  [8]  studied  longitudinal  waves  in  a  viscoelastic  rod 
caused  by  a  sinusoidal  stress  applied  at  one  end.  They  included  time  as 
an  independent  variable  and  obtained  time  histories  of  the  stress  and 
temperature  using  the  method  of  finite  differences.  Chang  [9]  has  stu¬ 
died  the  transient  temperature  profile  in  an  Infinite  viscoelastic 
medium  exterior  to  a  cylindrical  cavity  subjected  to  an  oscillatory, 
axially  symmetric  boundary  shear.  More  recently,  Knauss  [10],  Hege- 
mier  and  Morland  [11],  and  Edelstein  [12]  have  studied  thermal  effects 
in  viscoelastic  solids.  Knauss  considered  the  dynamic  response  of  a 
] sng  viscoelastic  bar  due  to  a  step  displacement  at  the  end;  Hegemier 
and  Morland  the  response  of  a  viscoelastic  half  space  to  the  sudden 
application  of  a  heat  source  distribution;  and  Edelstein  an  abalating 
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cylinder.  All  three  of  them,  however,  neglected  thermomechanical  coupl¬ 
ing,  i.e.,  the  heat  source  due  to  dissipation. 

Studies  of  thermal  effects  in  nonlinear  viscoelastic  media  have 
been  reported  in  [13]- [20].  In  [14]  Coleman  and  Gurtin  discussed  the 
velocity  and  the  growth  and  decay  of  one- dimensional  acceleration  waves 
in  conductors  and  in  media  which  do  not  conduct  heat.  In  [15]  some  of 
the  findings  of  [14]  ware  extended  to  three-dimensional  acceleration 
waves.  In  [16]  one- dimensional  shock  waves  in  non-conducting  viscoelas¬ 
tic  media  were  discussed.  In  [17]  Achenbach,  Vogel  and  Herrmann  dis¬ 
cussed  the  propagation  of  shock  waves  in  a  viscoelastic  material  with 
temperature  dependent  properties  and  thermomechanical  coupling.  The 
waves  were  caused  by  sudden  application  of  surface  tractions  and  the 
analysis  was  concerned  with  small  strains  and  changes  in  temperature. 
Oden  and  Ramirez  [18]  manipulated  their  relations  for  finite  difference 
application.  McGuirt  and  Lianis  [19]  compared  experimental  results  of 
finite  uniaxial  and  equal  biaxial  tests  for  non- isothermal  conductors 
with  theory  and  reported  good  agreement;  while  Lianis  [20]  suggested 
integral  constitutive  equations  for  nonlinear  thermoviscoelastic  mate¬ 
rials.  Random  temperature  effects  have  been  considered  by  Parkus  and 
Bargmann  [21]  m  a  recent  paper. 

This  report  is  concerned  with  efficiently  obtaining  steady  state 
solutions  for  stress  and  temperature  in  longitudinal  wave  propagation 
in  a  viscoelastic  rod  with  thermal  coupling  subjected  to  a  sinusoidal 
stress  at  one  end.  Huang  and  Lee  in  [8]  obtained  time  histories  for 
the  stress  and  temperature  for  this  problem.  They  set  up  the  govern¬ 
ing  nonlinear  partial  differential  equations  based  on  Schapery’s  model 
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and  solved  them  r  •  a  method  of  finite  differences.  In  practice,  the 
temperature  in  the  rod  increases  with  time  owing  to  a  continuous  dissi¬ 
pate  mi  of  mechanical  xr.ergy  into  heat  and  a  stf ady  state  might  be 
reached  if  the  rate  of  heat  gei  iration  by  the  dissipative  source  equals 
• '  raj  x  of  heat  f low  .ng  out  of  the  rod.  Such  a  steady  state  yields 
mos..  sev'  e  fmperature  conditions  vMch  are  of  primary  interest  in 
d  ->Ign.  Ir-  tu  .b  ases  it  is  ni  -h  more  efficient  to  obtain  the  steady 
state  direc..  .y  instead  of  following  the  complete  time  history  of  the 
process.  In  this  report  this  is  done,  by  solving  a  set  of  nonlinear 
ordinary  differential  equations  by  a  method  of  finite  differences.  A 
very  useful  algorithm  is  developed  by  'partial  linearization’  of  the 
equations  and  numerical  solutions  for  a  Lockheed  solid  propellant  are 
obtained  over  a  wide  range  of  driving  stresses  and  frecuencies.  For 
the  cases  studied,  steady  state  solutions  fox'  stress  and  temperature 
are  obtained.  The  very  high  temperatures  reached,  however,  would  cause 
failure  of  the  rod.  The  word  ’instability’  in  the  title  may  tht  bx 
interpreted  from  an  engineering  (practical)  point  of  view,  although  the 
existence  of  a  mathematical  instability  cannot  be  ruled  out  since  the 
complex  modulus  approaches  zero  as  the  temperature  increases.  In  fact, 
if  the  driven  end  is  insulated,  instead  of  being  kept  at  constant  tem¬ 
perature  as  in  our  problem,  unbounded  increase  of  the  compliance,  cor¬ 
responding  to  vanishing  modulus,  would  lead  to  unbounded  dissipation 
of  power  at  the  driven  end.  This  would  then  cause  infinite  tempera¬ 
tures  and  a  steady  state  temperature  solution  would  no  longer  exist. 
These  aspects  of  the  problem  are  discussed  in  more  detail  later  on  in 


the  report. 


FORMULATION  OF  THE  -PROBLEM 


The  formulation  of  the  p  Vbiem  closely  follows  that  in  [8]  with 
the  essential  difference  that  ’■no  steady  state  problem  is  considered 
here.  Variational  principles  were  used  by  the  author  to  obtain  a  nu¬ 
merical  solution  to  this  problem  for  one  value  of  driving  stress  and 
frequency  in  [25]. 

Let  us  define  stress,  strain  and  displacement  as  the  real  parts 


~  ^  imt  ~  iaat  ~  iu)t 

$  -  o  e  ,  e  =  ee  ,  u  =  ue 


where  i  <i)  is  the  frequency  and  t  is  time.  The  complex 

quantities  cr  %  e  and  u  ,  functions  of  the  space  variable  x  ,  will 
be  referred  to  simply  as  the  stress,  strain  and  temperature  respec¬ 
tively. 

Let  us  consider  a  viscoelastic  rod  of  length  &■  insulated  on  its 
lateral  surface,  as  in  Fig.  1.  The  left  end  is  free  while  the  right 
end  has  a  prescribed  stress  cos  ut>t  (cr^  real).  The  temperature  of 
the  vibrator  is  assumed  constant  at  Tq  while  a  radiation  boundary  con¬ 
dition  is  assumed  at  x  =  0  . 

The  boundary  conditions,  therefore,  are 


x  =  0 


c  =  0 


—  ~  c(T  -  T  ) 

u::  o 


T  =  T 


X  =  l 
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The  steady  state  energy  equation  in  one  dimension  for  the  cycle 
averaged  temperature  is  given  by  (see  [7,8,25]) 

d2T 

K  —5  +  2D  =  0  (7) 

dx 


Here  2D  is  the  cycle  averaged  value  of  the  mechanical  dissipa¬ 
tion  function 


iwl+°b 


where  Re  denotes  the  real  part  of  the  complex  argument. 

We  note  here  that  part  of  the  mechanical  dissipation  is  trans¬ 
formed  into  heat  while  the  remaining  causes  increase  of  potential  energy 
Heat  is  also  produced  due  to  dilatational  compression.  These  coupling 
terms  due  to  dilatation  and  potential  energy,  however,  are  periodic 
and  drop  out  if  the  energy  equation  is  integrated  over  a  cycle.  Thus 
the  temperature  T  here  is  the  cycle  averaged  value  of  the  true  tem¬ 
perature.  T  is  assumed  to  attain  a  steady  state  after  a  sufficiently 
long  time  and  become  independent  of  time.  Strictly  speaking,  the  tem¬ 
perature  has  small  cyclic  variations  about  a  mean  value  as  a  result  of 
the  cyclic  variations  of  the  potential  energy,  dilatation  and  dissipa¬ 
tion,  but  these  are  quite  negligible. 

We  note  also  that  while  the  mechanical  transients  vanish  within 
one  or  two  cycles  as  reported  in  [A],  thermal  transients,  i.e.,  the  slow 
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increase  of  temperature  with  time,  can  continue  over  several  hours  be¬ 
fore  it  settles  down  to  a  steady  state,  if,  in  fact,  a  steady  state 
exists.  In  the  numerical  examples  given  in  this  report,  such  a  steady 
state  is  found  to  exist  and  is  obtained  directly  without  integrating 
forward  in  time. 

Equations  (5),  (6),  and  (7)  constitute  a  set  of  three  nonlinear 

ordinary  differential  equations  for  the  unknowns  ,  a and  T  . 

The  associated  two  point  boundary  conditions  are  given  by  Eq.  (2). 

Thermomechanical  coupling  is  caused  by  the  nonlinear  dissipative 

source  D  .  The  real  and  imaginary  parts  of  the  complex  compliance 
★ 

D  are,  in  general,  nonlinear  functions  of  temperature. 

For  a  Lockheed  solid  propellant  over  a  wide  reduced  frequency 
range  we  can  express  and  by  the  following  empirical  formulae 

(see  [8]) 


B  v 

DL  =  ctu)  (T  -  Tx)  Y 

D2  =  c2uT(T  -  Tt)V 


(8) 


where  ,  c2  *  P  ,  Y  and  T^  are  constants. 

Let  us  define  the  nondimens ional  quantities 

x  T  -  T  _  i/j 

q  -  1  »  T  =  -T  "-_~t  *  *  =  ,  X  =  [2Kojp(T  -  T^]  *  , 

o  1 


s  ~  xa  ,  so  =  XaQ  ,  Sl  =  Xc^  ,  s2  =  Xa2 


d^  =  £2ai2pD1  ,  d2  =  aj2pD2 


(9) 
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The  governing  differential  equations  (5) ,  (6) ,  and  (7)  reduce  to 
the  simplified  form 


s'^  +  +  djS^  =  0 


s2  +  dlS2  '  d2Sl  “  0 


T"  +  d2(s^  +  s2)  =  0 


(10) 

(ID 

(12) 


where 


and 


v  Y 

-  c^aT  ,  T 


a  «=  u)^+^(Tq  -  T^)\  a  constant 


(  )*  s 


dq 


The  boundary  conditions,  Eq.  (2),  become 


q  =  0  :  s1  =  0  ,  s0  =  0  ,  t'  =  k(t  -  1) 


q=l:s1=sQ  ,  s2=0  ,  t=1 


(13) 


We  thus  have  a  two  point  boundary  value  problem  which  is  solved 
by  a  method  of  finite  differences  (see  [22,23,24]). 


METHOD  OF  SOLUTION 

1.  Partial  Linearization 

Equations  (10)  and  (11)  are  linear  in  the  stresses  given  the  tem¬ 
perature  distribution.  However,  Eq.  (12)  is  not  linear  in  temperature, 
given  the  stresses,  since  d2  is  a  nonlinear  function  of  temperature. 
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If  Eq.  (10)  Is  multiplied  by  s2  ,  Eq.  (11)  by  ,  and  the  resulting 
second  equation  subs trac ted  from  the  first,  we  get 

s2s’^  -  +  d2(sj  +  s2)  =  0  (14) 

and  using  Eqs.  (14)  and  (12)  we  have 


T"  =  S„S 


2  1 


S,  s 


12 


(15) 


Integrating  once  with  respect  to  the  independent  variable  q 


T'  -  t’(0)  =  s2s|  -  s1s^ 


(16) 


Now  Eq.  (16)  is  linear  in  temperature  given  the  stress  distribution 
and  can  be  used  to  replace  Eq.  (12)  in  the  original  set  of  Eqs.  (10), 
(11)  and  (12).  Thus,  the  new  set  of  governing  differential  equations 
((10),  (11),  and  (16))  are  'partially  linear'  and  this  set,  together 
with  the  associated  boundary  conditions  (Eq.  (13))  is  used  to  obtain 
solutions  for  the  stress,  strain  and  temperature  distributions. 

2.  Finite  Difference  Algorithm 

The  iterative  procedure  works  as  follows: 

Guess 


,  Equations  of  Motion 

=>  (10)  and  (11) 

+  boundary  conditions 


,(j> 


,<j) 


Modified  Energy 
=>  Equation  (16) 

+  boundary  conditions 


tU+D 


where  the  superscript  j  denotes  the  jth  approximation  to  the  solution. 


Iteration  is  continued  till  convergence  is  achieved. 


Let 
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The  boundary  conditions  give 


0 

s 

11 

»—< 

0 

0 

0 

Let  V.  be  a  sequence  of  2x2  transfer  matrices  defined  as 


\*VW  CK-0,1,...,N) 


I 

V  = 
o 


0  0 


0  0 


and  substituting  Eq.  (19)  into  Eq.  (17)  we  get  the  recurrence  relation 

.-  -r  1 

vk  =  -  L\  +  ViJ  (k  =  1,2, . . .  ,N)  (21) 

Now  with  V  known  from  Eq.  (20),  Eq.  (21)  gives  V,  ,  V.  ,  ..., 
o  i  £■ 

VN  and  then  Eq.  (19)  gives  »  t'n.1  »  •  •  •  >  T'1  • 

Thus,  the  equations  of  motion  are  completely  solved  for  a  certain 
guess  of  the  temperature  distribution.  The  use  of  transfer  matrices 
is  a  big  advantage  since  only  2x2  matrices  need  to  be  inverted. 

The  modified  energy  equation  (16)  and  its  associated  boundary  con¬ 
ditions  become 


<sm  -  s'-i>  sl 


^Sk+l  "  Sk-l'  Sk 


(Tk+l  "  Tk-1) 


-  H(To  -  1)  (k  =  1,2,... ,N) 


•T„  +  4T.  -  3T  =26  K(T  -  1) 
2  1  o  o 


Vi  =  1 


Tk  =  ^kTo  +  Yk  ^  =  0,1 . IW'1) 


“k-sLi-sli  <k  =  i-2 . N) 


sk  -  sL  -  sL  <k  -  *'2 . » 


Substitution  of  Eq.  (23)  into  Eq.  (22)  leads  to 


B  =  1  ,  Bj_  =  1  +  6k 


- 0  •  \  -  \  <SX  -  s?zi>  -  6k 


and  the  recurrence  relations 


Sk+1  =  \-l  +  26K  (k  -  1,2,...,N) 


vx  -  Vi  -  26r  +  vi  -  \sk 


and  once  &k  »  \  (k  ®  0,1, . . . ,N+1)  are  known,  we  have 


T  =  — ->j2  t  =  ! 

<Vi  ’  ** 


and  (k  =  1,2, . . . ,N)  are  obtained  from  Eq.  (23). 

Thus,  the  modified  energy  equation  (16)  is  solved  for  the  tempera 
ture,  given  the  stresses. 


The  iterative  procedure  outlined  in  the  beginning  of  this  section 
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is  allowed  to  proceed  till  successive  approximations  for  are  within 

certain  prescribed  convergence  limits.  The  accuracy  of  the  solution  is 
verified  by  back  substitution  into  the  original  finite  difference  equa¬ 
tions.  Different  mesh  sizes  are  chosen  in  such  a  way  that  halving  the 
step  size  causes  a  change  of  .  57o  or  less  in  the  results. 

Thus,  the  steady  state  values  of  the  real  and  imaginary  parts  of 
the  nondimens ional  stress,  s^  and  s2  and  the  nondimens ional  tempera¬ 
ture  t  are  obtained.  The  nondimens ional  stress  at  any  time  is 

Re(s  eltUt)  =  s^  cos  tut  -  s^  sin  tut 


RESULTS  AND  CONCLUSIONS 

Numerical  calculations  have  been  carried  out  for  the  following 
data  for  a  Lockheed  solid  propellant  [8] 

cx  =  4.61  x  10"11  (psi)_1(sec)VF}‘Y 

c2  =  1.62  x  10"11  <psi)"1(sec)P(0F)"Y 

3  «  -  0.214  ,  y  =  3.21 

H  =  1.0  ,  T  =  65°F 

o 

Tl  =  -125°F  ,  l  =  3  in. 

42p  =  1.023  x  lo“4  pst-sec2 

2Kp(io  -  tl)  =  8.08  x  10-4  psi2-sec 
In  [8],  K  should  read  1.0  instead  of  0.1. 


Various  values  for  the  frequency  u)  and  the  driving  stress  sq 
were  used. 

Figure  2  shows  the  value  of  t  as  a  function  of  q  for  various 

values  of  the  driving  stress  s  .  The  frequency  u)  is  kept  constant 
4 

at  10  radians /second.  The  steady  state  temperature  profiles  are 
found  to  exist.  However,  the  temperature  response  is  markedly  non¬ 
linear,  increasing  much  more  rapidly  than  stress  and  reaching  danger¬ 
ously  high  levels  for  low  values  of  driving  stress.  The  constitutive 

*  * 

law  used  here  assumed  that  the  complex  Young's  modulus  E  =  1/D  ap¬ 
proaches  zero  as  the  temperature  T  -*•  “  (see  Eq.  (8)).  Schapery  [6] 

has  an  experimental  graph  for  a^,(T)  for  our  Lockheed  propellant  which 

* 

shows  that  for  large  temperatures  E  0  for  a  finite  value  of  T  . 

* 

E  — 5>  0  means  zero  wave  velocity  in  the  material  and  consequently 
failure.  Thus,  the  temperature  levels  attained  in  our  calculations 
(as  shown  in  Fig.  2)  would  lead  to  melting  and  failure  of  the  rod,  even 
though  a  mathematical  steady  state  solution  exists  and  T  is  finite 
everywhere . 

Other  interesting  observations  can  be  made  regarding  Fig.  2.  The 
rapid  softening  of  the  material  of  the  rod  leads  to  large  values  of 
strain  (around  20%  near  the  driven  end  for  sq  =  3.0)  so  that  further 

calculations  for  larger  values  of  s  would  not  be  valid  within  the 

o 

realm  of  linear  viscoelasticity.  With  increasing  stress  the  location 
of  the  temperature  maximum  asymptotically  approaches  the  driven  end  of 
the  rod  since  the  larger  damping  at  higher  temperatures  leads  to  rapid 
attenuation  of  the  stress  away  from  the  driven  end.  Consequently,  most 
of  the  mechanical  energy  is  dissipated  as  heat  near  the  driven  end  of 


Pi 
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the  rod.  Another  consequence  of  softening  of  the  rod  is  lower  stress 
wave  velocities  and  therefore  smaller  wave  lengths  since  the  driving 
frequency  is  kept  the  same.  These  effects  on  the  s^  and  profiles 

are  clearly  evident  in  Figs.  3,  4  and  5  which  show  s^  and  s£  as 

4 

functions  of  q  for  s  =  0.5,  1.5,  and  2.5  with  U)  =  10  radians/ 

o 

second.  In  Fig.  5  the  stress  is  practically  zero  beyond  about  a  third 
of  the  length  of  the  rod  from  the  driven  end. 

The  effect  of  driving  frequency  on  the  maximum  temperature  is 
shown  in  Fig.  6  for  two  values  of  the  driving  stress.  The  maximum 
temperature  is  seen  to  have  several  peaks  at  certain  critical  frequen¬ 
cies.  Increase  of  driving  stress  gives  higher  and  more  frequent  peaks 
at  different  critical  frequencies.  This  is  because  more  input  of 

mechanical  energy  leads  to  more  dissipative  heating  and  therefore 

* 

larger  temperatures.  This,  in  turn,  lowers  E  and  the  change  of 
properties  of  the  rod  changes  the  critical  frequencies.  Also,  softening 
of  the  material  lowers  the  stress  wave  velocity  in  the  rod  and  this 
leads  to  a  lower  fundamental  frequency  for  higher  driving  stress. 

The  effect  of  driving  stress  is,  of  course,  largest  near  these 
critical  frequencies.  As  an  example,  the  temperature  and  stress  plots 
for  (U  =  2500  radians /second  (the  fundamental  critical  frequency  for 
a  =2.84  psi)  are  shown  for  o  =1.42  psi  and  a  =2.84  psl  in 
Figs.  7,  8,  and  9.  A  very  large  temperature  rise  is  obtained  for  a 
doubling  of  the  driving  stress  and  the  and  S£  plots  are  consid¬ 

erably  different  for  «  2.84  psi  compared  to  those  at  1.42  psi.  Thus, 
these  critical  frequencies  are  of  utmost  importance  in  design. 

An  attempt  was  made  to  compare  the  critical  frequencies  of  the 
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coupled  problem  (as  in  Fig.  6)  with  the  natural  frequencies  of  the  un¬ 
coupled  problem  where  the  temperature  distribution  in  the  rod  is  al¬ 
ready  specified.  If  a  certain  temperature  distribution  is  prescribed 
in  the  rod,  we  have  an  inhomogenious  viscoelastic  rod  where  the  complex 
Young's  modulus  is  known  as  a  function  of  the  distance  along  it.  We 
thus  have  a  linear  problem  and  the  complex  natural  frequencies  of  the 
rod  can  be  determined  by  solving  the  resultant  eigenvalue  problem. 

This  was  done  using  a  method  of  finite  differences.  The  details  of  the 
method  are  given  in  the  appendix. 

Let  the  temperature  distribution  from  the  coupled  problem  for 
oq  =  2.84  psi  ,  u)  *  2500  radians /second  be  prescribed  in  the  rod. 

The  table  gives  a  comparison  of  the  first  three  uncoupled  natural  fre¬ 
quencies  with  the  critical  frequencies  of  the  coupled  problem  from 
Fig.  6. 


Uncoupled  Frequency 
Radians / second 

Coupled  Frequency 
Radians / second 

“l 

1739.04  +  333.031 

2500 

U) 

2 

4195.78  +  803.501 

5000 

u> 

3 

6708.83  +  1284.751 

7500 

Also,  the  fundamental  uncoupled  frequency  in  case  the  temperature 
distribution  is  prescribed  for  =  1.42  psi  ,  u>  *  4000  radians/ 
second  is  3484.36  +  667. 26i  radians /second  compared  to  the  coupled 
fundamental  critica'  •  -equency  of  4000  radians /second. 

The  response  of  the  uncoupled  rod  to  forced  oscillation  at  various 
frequencies  is  shown  in  Fig.  10.  Here  the  prescribed  temperature 
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distribution  is  that  obtained  from  the  general  coupled  problem  for  gq 
2.84  psi  ,  u)  =  2500  radians /second.  We  see  from  Fig.  10  that  owing  to 
the  presence  of  damping  the  peak  is  shallow  and  occurs  around  1800 
radians /second  (compare  with  <u^  in  table).  Thus  the  eigenvalues  of 
the  uncoupled  problem  give  us  an  idea  of  those  of  the  coupled  problem 
and  shed  some  light  on  the  correlation  between  mechanical  and  thermal 
resonance . 

This  study  shows  that  mechanical  failure  of  viscoelastic  polymers 
under  cyclic  loading  can  occur  at  very  low  values  of  driving  stress, 
especially  if  the  driving  freqeuncy  is  close  to  one  of  the  critical  fre 
quencies  of  the  coupled  problem.  For  purposes  of  design,  the  direct 
determination  of  the  steady  state  solutions  for  temperature  and  stress 
appear  superior  to  the  method  used  by  Huang  and  Lee  in  [8]  where  the 
complete  time  histories  of  the  above  mentioned  quantities  were  deter¬ 
mined.  A  lot  of  computational  effort  is  saved  since  the  method  used 
here  requires  solution  of  ordinary  rather  than  partial  differential 
equations  in  a  computer.  'Partial  linearization'  of  the  problem  leads 
to  a  very  useful  algorithm  which  converges  rapidly  for  a  wide  range  of 
driving  stresses  and  frequencies. 


-19- 


APPENDIX 

Response  of  a  Linear  Inhomogenious  Viscoelastic  Rod  with  the 
Temperature  Distribution  T(q)  Prescribed 

1.  Determination  of  Natural  Frequencies 

The  equations  of  motion  of  the  rod  can  be  written  as 

s"  +  ds  *  0  (A-l) 


where 


s  *  Sj^  +  is2  ,  d  *=  d^  -  id2 

d(q)  «  (cx  -  ic2)(To  -  T^y  £2p(T(q))Y  u>2+B 


Let 


d(q)  =  r(q)  p 


where 


P  -  ID 


2+3 


p(q)  =  l/r(q)  “ 


(cL  +  ic2) 


(C1  +  c2)  (To  '  T1)V  £2P<t^»Y 


Now  Eq.  (A-l)  becomes 


p(q)  s'*  +  ps  *  0  (A- 2) 

s(0)  »  s(l)  »  0 


which  is  an  eigenvalue  problem  with  p  ,  p  and  s  complex. 
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,  <lk  =  k6  (k  =  0,1 . H+l) 

s<1k>  =  sk 

p(qk)  =  Pk  (k  =  0,1,... ,Nfl) 


we  have,  from  Eq.  (A- 2) 

-Pk  +  !Vk  '  pk5k+i  ■ p52  sk  (k  • l-2 . N) 

s0  =  sm  =  0 

In  matrix  form,  this  can  be  written  as 

Vs_  =  n&2  £ 

2p^  -p^  0  0  ....  0  0  0 

-p2  2p2  -p2  0  -  0  0  0 

0  _p3  2p3  ~P3  -  0  0  0 

0  0  0  0  -  ”PN-1  2pN-l  "PN-1 

0  0  0  0  ....  0  ”PN  2pN 
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Thus,  v.  (i  =  1,2, 
natural  frequencies 


. ,N)  are  the  eigenvalues  of  P  ,  the 


<o. 

l 


,v 


-ti 


1 

2+5 


(i  =  1,2,. ...N) 


We  know,  from  [22]  that  these  numbers  approximate  the  N  smallest 
natural  frequencies.  Also,  if  arranged  in  an  increasing  order  of  magni¬ 
tude,  the  approximation  is  very  good  foi  the  first  natural  frequency 
but  the  accuracy  of  the  approximation  decreases  rather  rapidly  with  in- 
creas ing  N  . 

In  order  to  determine  the  complex  eigenvalues  of  the  complex  ma¬ 
trix  P  ,  it  is  enough  to  solve  for  the  eigenvalues  of  the  real  matrix 
P^  ,  since 

p  -  h  +  ip2  -  pi  +  ^r) 


and,  of  course,  the  eigenvalues  of  P  , 

ic_ 

v .  =  ( 1  + - )  a . 

i  '  c^f  i 

where  are  the  eigenvalues  of  P^  . 


Since  Re(r(q))  >  0  for  all  q  ,  all  are  real  and  positive 

(see  [22]). 

The  normal  modes  can  be  determined,  if  desired,  by  finding  the 
eigenvalues  of  P  . 


2-  Response  to  Forced  Oscillations 

This  problem  is  a  special  case  of  the  general  problem  with  thermo¬ 
mechanical  coupling  and  can  be  easily  solved  using  a  slight  modification 
of  the  algorithm  for  the  coupled  problem. 
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Figure  1.  Boundary  conditions  for  the  problem. 
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